Quantum percolation and Anderson transition point for transport 

of a two-state particle 
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Quantum percolation describes the problem of a quantum particle moving through a randomly 
frozen medium. While certain similarities to classical percolation exist, the dynamics of quantum 
percolation has additional complexity due to the possibility of Anderson localization. Here we 
show that this strongly influences the percolations threshold by considering a directed two-state 
quantum walk on a two dimensional space. To do this we determine the Anderson transition 
point (the quantum equivalent to the classical percolation threshold) for three fundamental lattice 
geometries (finite square lattice, honeycomb lattice, and nanotube structure) and show that it 
differs significantly from the classical value and tends towards unity for increasing lattice sizes. 
Beyond the fundamental interest for understanding the dynamics of a two-state particle on the 
lattice (network) with disconnected vertices, our study also sheds light on the transport dynamics 
in various quantum condensed matter systems and the construction of several quantum information 
processing and communication protocols. 
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Percolation theory is a mathematical model originally 
developed to describe the dynamics of particles in ran- 
dom media [1]. Due to its fundamental nature it has since 
established itself as an area of research on its own right 
and found numerous applications in diverse fields. These 
include fluid dynamics, fire propagation, many body sys- 
tem in classical and quantum physics, information theory, 
dynamics in biological system and in chemical compo- 
nents [2, 3]. 

The main figure of merit which determines the trans- 
port efficiency of a particle in percolation theory is the so- 
called percolation threshold [4] . To illustrate its meaning 
in the classical setting, one can consider transport on a 
square lattice with neighbouring vertices connected with 
probability p. When p = 0, all vertices are disconnected 
from each other and no path for the particle to move 
across the lattice exists. With increasing p more and 
more vertices will be connected and when p = p c = 0.5 
the first connection across the full lattice is established. 

The corresponding problem of percolation of a quan- 
tum particle differs from the classical setting in that it 
contains an additional degree of freedom. In a disordered 
system the interference of the different phases accumu- 
lated along different routes during the evolution can lead 
to the particle's wave function becoming exponentially lo- 
calized. This process known as Anderson localisation [5- 
7] and has recently been observed in different disordered 
systems [8-10]. Quantum interference therefore becomes 
as important in quantum percolation as the existence of 
the connection between the vertices, making it a more 
intriguing setting when compared to the classical perco- 
lation [11-14]. 

Transport of a two-state quantum system across a large 
network is an important process in quantum information 
processing and communication protocols [15] and by to- 
day many physical systems are tested for their scalability 
and engineering properties. Furthermore, in last couple 
of years quantum transport models have also shown a cer- 



tain applicability to understanding transport processes in 
biological and chemical systems [16, 17]. Since these nat- 
ural or synthetic systems are not guaranteed to have a 
perfectly connected lattice structure, it is important to 
consider the possible role quantum percolation can play 
in understanding transport in these systems. To model 
the dynamics of the quantum particle we choose the pro- 
cess of quantum walks, which in recent years have been 
shown to be an important and highly applicable mech- 
anism [18]. Recently, first studies on quantum walks in 
percolating graphs have been reported for circular and 
linear geometries [19] and for square lattices using a four- 
state particle [20] . 

Here we present a new model of a directed discrete- 
time quantum walk (D-DQW) to study quantum percola- 
tion of a two-state particle on a two-dimensional (2D) lat- 
tices. The D-DQW is a physically relevant model which, 
for example, describes the dynamics of transverse propa- 
gation of electrons and photons in 2D lattices and optical 
waveguides. Our main finding is that the possible locali- 
sation due to quantum interference at disconnected ver- 
tices leads to an Anderson transition point (p a ) which dif- 
fers significantly from its classical counterpart. To show 
this we consider three lattice geometries (square lattice, 
honeycomb lattice and nanotube geometry) and numeri- 
cally determine p a for different lattice sizes. We also show 
that generalising the D-DQW dynamics to allow for more 
degrees of freedom results in the same values for p a , hint- 
ing at the significance of D-DQW as a fundamental and 
efficient model to describe dynamics. Interestingly, we 
find that p a is lattice size dependent and tends towards 
unity for larger lattices. This is in contrast to the classi- 
cal case, where the percolation threshold is independent 
of the lattice size. We also find that p a is smaller for fi- 
nite size honeycomb lattices and the nanotube geometry 
when compared to the square lattice. 

This model can effectively be used as a framework to 
explore the dynamics for different evolution protocols and 
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find the most suitable model to describe the quantum 
transport in various natural physical systems. For engi- 
neered systems one can explore the dynamical processes 
for given lattice structures of different size which might 
have a large degree of disconnected vertices and test their 
suitability for quantum information storage or quantum 
communication protocols. 

Model - Let us first define the dynamics of a D-DQW 
on a completely connected square lattice. The Hilbert 
space of the particle, Hp, is represented by its internal 



states, | 1) 



and | |) = 



and the Hilbert space 



of the square lattice of dimension n x n, Wl, is repre- 
sented by the position of the vertices, \x,y). Each step 
of the walk consists of first evolving the particle into a 
superposition of the internal states using a coin flip op- 
eration We followed by the shift in the x and y directions 
using W x and W£, respectively. The index c stands for 
completely connected vertices in the lattice. The Wg and 
W£ are the same operations as used in the standard ID 
DQW [21-23], 
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For the shift in y-direction, we define the position de- 
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where a x _ y = a x +i,y and f3 x _ y = P x -x,y to ensure trans- 
port in the positive y-direction only. The state of the 
particle after t steps is then given by 



(4) 



x y 



and we assume the initial state to be \*f?i n ) = 
(cos((5/2)U) + e"'sin((5/2)| t)) ® ||n/2"|,l) [24] and 

\lpx,y) = (ax, y \l)+0x,y\t})®\x,y) = Vi„ 

In Fig. 1(a), we show the path taken by a particle on 
a perfectly connected square lattice. The probability of 
detecting the particle outside of the lattice after t steps 
is given by, 



P(t) = l-^2^2(x t y\p(t)\x,y), 



(5) 
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where p(t) = |4 r t )(^ r t |. Note that due to the directed 
component of the walk one finds P(t) = 1 for t > n. 
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FIG. 1: Schematic of the path taken by a two-state particle 
on a square lattice. Green, red and blue arrows represent the 
direction of | J,), f) and both the states, respectively, when (a) 
all vertices are perfectly connected and (b) some connections are 
missing leading to localization at the highlighted positions. 



An alternative description of this dynamics can be 
found in the form of a differential equation by consider- 
ing the state of the particle at vertex (x, y) as a function 
of 6 at any time t, 



^ = 008(0)^+^.! - i M^l-l,y-l (6) 

4>l v = cosC^vLi,^! - i MWi+i, v -i ■ (7) 

These equations can be easily decoupled and written as, 
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Subtracting then 2 [1 + cos(0)J ipx)y from both sides of 
Eq. (8), we obtain a difference form which can be written 
as a second order differential wave equation 
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One can see that the wave-packet propagates along both 
directions on the x-axis at a rate of D x = \ ^ cos(6)\ per 
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FIG. 2: Increase in the percolation probability and the 
Anderson transition point as a function of the percentage 
of connections, the lattice size and the coin parameter 9. 

Probability for a particle to be found outside of the square lattice 
as a function of the percentage of connected vertices for different 
lattice sizes for 9 = n/A. The Anderson transition point can be 
seen to approach unity for increasing lattice size and also has a 
strong dependence on 9 for 9 < n/A (see inset). 



time step [25] and in the positive y direction at the rate 
of D v = 1. For a completely connected lattice the total 
displacement after time t in the x- and the y-direction 
will therefore be D x (t) = \t<J caa(ff)\ and D y (t) = t, re- 
spectively. 

Let us now consider lattice structures in which some of 
the edges connecting the vertices are missing. If all edges 
connecting the vertices on the lattice can be represented 
by a set E, the shift operator along the x- and y-axis can 
be written as 



W x 



if ({x,y},{x + l,y} 
and ({x,y},{x- l,y} 



G E 
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where W° and Wy are given in Eqs. (2) and (3), and 

W A X = £ [(I I) (I | + | t)(t I) ® \x,y)(x,y\] (12) 
x,y 

W v = J2 [(|-*,v><-*,vl + l+*,»>(+»,vl) ® \x,v)(x,v\ ■ 

(13) 



Here the index d stand for vertices with disconnected 
edges. In Fig. 1(b), we show an example of a path taken 
by a particle in a lattice where some connections are miss- 
ing. In a classical setting the percolation threshold for 
this square lattice can be calculated to be p c = 0.5 and is 
known to be independent of the lattice size. In a quan- 
tum system, however, a disconnected vertex breaks the 
well defined interference of the multiple traversing paths 
and results in trapping a fraction of the amplitude at 
the disconnected vertex. This mix up of the interference 



due to the disorder is known to result in Anderson lo- 
calization [5, 6] and consequently a large percentage of 
connected vertices are required to reach a non-zero proba- 
bility for the particle to cross the lattice. In the following 
we will call the probability for the particle to cross the 
lattice the percolation probability, CCf); an d its critical 
value at which the percentage of connected vertices p is 
large enough to reach £(p) = 0.01 the Anderson transi- 
tion point, p a . It is obtained numerically by averaging 
over several runs and one of our main findings is that it 
is significantly larger than p c for the same lattice. We 
also find that it converges towards unity with increasing 
lattice size. To show this we plot the average £(p) as a 
function of p for square lattices of different size in Fig. 2. 

When all vertices in the square lattice arc fully con- 
nected, initial position (x, y) = (\n/2\ , 1) and 9 = 0, the 
two basis states move away from each other along the x- 
axis (no interference takes place) and exit from the sides 
when half way through the y-axis. For finite values of 9, 
the exit point is pushed towards the positive y direction 
due to interference in both directions along the x-axis. In 
the absence of perfect connections between the vertices 
in a lattice of finite size, the paths are altered resulting in 
interference for all values of 9. For small values of 8 how- 
ever, a large fraction of the states still exits along the 
sides without interference and therefore only a smaller 
number of connections play a role. This explains the in- 
crease in p a as a function of 9 and to show this we plot 
p a as function of 9 in the inset of Fig. 2. For large lattice 
size, dependence of p a on 9 will be negligible. 

In the same way as for Eq. (9), a general differential 
equation for one step evolution on a lattice with discon- 
nected vertices can be written as [26], 



p 3 + 2p\l-p)+p{l-pf 
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+ [1 - cos(f5)] \2p 3 + 5p 2 (l -p) + 8p(l - pf + 2(1 - pf 



-p 2 (l -p)[l - 3cos(0)] 



V4, = o. 



(14) 



When p = 1 the expression reduces to Eq. (9) and when 
p = all the derivative terms vanish. Here the terms con- 
taining the first and second order derivatives contribute 
to the propagation and the other terms describe local- 
ization. Even for values of p — > 1 one can see the dom- 
inance of the non-derivative terms, which explains the 
large value for p a . 

Generic extension: The assumption of having the same 
coin operator at each lattice site is a rather strong one 
and in the following we will relax this condition to ac- 
count for applications in more realistic situations. For 
this we replace 9 by a vertex dependent parameter, 
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FIG. 3: Smaller spread in probability distribution for the 

D-DQW for 9 x , y when compared with 6 = 7r/4. Spread in the 
probability distribution in x-direction during propagation in the in- 
direction for a particle walking on a completely connected square 
lattice of size 100 X 100. In (a) 8 = 7r/4 is used for each step of the 
D-DQW and in (b) the coin is position dependent, 9 x ^ y S [0,tt]. 
The smaller spread in x direction for 6 x ^ y is clearly visible (red). 



FIG. 4: Increase in percolation probability and Anderson 
transition point as a function of the percentage of connec- 
tions and the lattice size. Percolation probability as a function 
of the percentage of connected vertices for square lattices of differ- 
ent size. Each vertex has a different value of 9 x ^ y , which has been 
randomly chosen from [0, it]. The p a shifts towards unity with in- 
creasing lattice size (see inset). 



@x,y G [Oj 71 "] an d note that ca&(6 x ,y) will be negative 
for any 9 XjV G [tt/2,vt]. This can be interpreted as the 
displacement of the left moving component to the right 
and the right moving component to the left along the x- 
direction, which in turn can lead to localization in trans- 
verse direction [27]. We show this in Fig. 3, where a typ- 
ical probability distribution of the particle in x-direction 
during its propagation in y-direction on a completely con- 
nected square lattice for a single realization. The average 
percolation probability as a function of p for this evolu- 
tion is then shown in Fig. 4 and, interestingly, we find 
that the disorder in the form 9 X , V does not result in any 
noticeable change in the value of p a when compared to 
6 = tt/4. This is due to nearly the same degree of in- 
terference [27] for both 9 XtV and 9 — ir/4 and highlights 
the dominance of the localization effects on disconnect 
vertices. This can also be seen from Eq. (14), where re- 
placing 9 with 9 x y will not result in significant change in 
the dynamics due to the dominance of p for p < 1 . 

Quantum percolation on a honeycomb lattice and 
nanotubes- Transport processes on honeycomb lattices 
and nanotubes have attracted considerable attention in 
recent years [28] and the two-state quantum percolation 
model can be expected to give useful insight into the be- 
haviour of quantum currents and their transition points. 
In Fig. 5, we show the path taken by the two-state D- 
DQW on the honeycomb lattice of dimension 9x9. Each 
step of D-DQW now comprises of the sequence W y W q We , 
where the quantum coin operator We and the operator 
for the directed evolution in the y-direction, W y , are same 
ones as used for the evolution on square lattice. Due to 
the honeycomb geometry, the transition W q from q to 
q ± 1 corresponds to a shift along two edges, first in the 
±x direction and then along the positive y direction. In 
Fig. 6(a) we show the average percolation probability as 
a function of the percentage of connected vertices with 
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FIG. 5: Schematic of the path taken by a two state particle 
on a honeycomb lattice. D-DQW on an honeycomb lattice of 
dimension 9x9, where (a) shows the situation where all connections 
between the vertices are present and (b) the situation where a finite 
number of connections are missing. This leads to localization at the 
marked position. Red and green represents the path taken by the 
two internal states and blue represent the path taken by both the 
states. 



randomly assigned value of 9 XyV £ [0, tt]. Similarly to 
the square lattice we find that the Anderson transition 
point p a is significantly larger than the classical perco- 
lation threshold p c = 0.652 [29] and also lattice size de- 
pendent. Note that compared to a square lattice of the 
same size, p a for honeycomb lattice is smaller, which gives 
the honeycomb structure a edge over the square lattice 
for quantum percolation using D-DQW. This is the re- 
versed situation compared to the classical scenario. This 
can be understood by considering the geometry of the 
honeycomb lattice: the particle needs to only propagate 
along two edges to shift one position in x direction and 
two positions in y direction, whereas in the square lattice 
transport along three edges is necessary to achieve this. 
A natural extension of the honeycomb lattice is to 
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TABLE I: Anderson transition point for finite systems. 
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FIG. 6: Increase in percolation probability and Anderson 
transition point for a honeycomb lattice and a nanotube 
geometry. Percolation probability as a function of percentage of 
connected vertices for (a) honeycomb lattices and (b) nanotube 
structures of different sizes. For the particle transport process the 
value of 9 has been randomly assigned from [0, it] for each vertex. 
With increase in lattice size, p a shifts towards unity. For nanotubcs 
of size n X y the transition point can be seen to be independent of 
n 



consider periodic boundary conditions in the x-direction, 
which transforms the honeycomb lattice into a nanotube 
geometry. This corresponds to transitions from q to 
(<Z±1 mod n), where n is the number of vertices along the 
x-axis and in Fig. 6(b) we show the percolation proba- 
bility for such a structure as a function of the percentage 
of connected vertices with randomly assigned values of 
x , y 6 [0,tt] at each vertex. One can see that the Ander- 
son transition point is same as that to the flat honeycomb 
structure with same number of edges in the transverse di- 
rection, which can be understood by realising that the pe- 
riodic boundary conditions make the particle encounter 
the the disconnected vertices more than once. The nan- 
otube with small number of vertices in radial direction 
therefore corresponds to an effectively larger flat system 
with the same defect density and from the earlier studies 
we know that larger lattices have higher Anderson tran- 
sition points. Due to the absence of the exit point along 
the radial axis, the only direction the particle can exit is 
the positive y-direction, which explains the independence 
of p a from the number of radial vertices. To summarise 
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our main results, we show a comparison of p a for the 
different geometries discussed above in Table I. 

In conclusion, we have investigated quantum percola- 
tion using a directed two-state DQW to model the quan- 
tum transport process. We have shown that the transi- 
tion point p a , beyond which quantum transport can be 
seen is relatively large when compared to the classical 
percolation threshold, p c . In addition, for finite lattice 
size and unlike the classical case, we have found that p a 
is size dependent and tends towards unity with increase 
in size. This suggests that even a small disconnected 
vertices in the large system obstructs the quantum trans- 
port. 

Comparing different lattice geometries we found that 
p a is smaller for honeycomb structures and nanotube ge- 
ometry as compared to square lattice. This variation 
suggests that one can explore the dynamics on different 
lattice structures to find one most suitable for a required 
purpose. For example, higher p a can allow for storage 
and lower ones for more efficient transport processes. 

The two-state quantum percolation model using D- 
DQW for transport process is a realistic model that can 
be used to study transport process in various directed 
physical systems such as photon dynamics in waveguides 
with disconnected paths or quantum currents on nan- 
otubcs. We have demonstrated its generality by allowing 
the parameter 9 to vary randomly at each vertex and 
showed that this does not lead to any significant change 
in p a . Given the current experimental interest and ad- 
vances in implementing quantum walks in various phys- 
ical systems [30] , we believe that our discrete model is a 
strong candidate for upcoming experimental studies and 
its differential equation form we have presented will also 
be interest for further theoretical analysis. 
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I 

SUPPLEMENTARY MATERIAL : QUANTUM PERCOLATION AND ANDERSON TRANSITION 
POINT FOR TRANSPORT OF A TWO-STATE PARTICLE 

For a completely connected lattice, it is straight forward to write the states of the particle at position (x, y) as 
function of 9 at any time t (see Eqs. (6) and (7) in the paper) and obtain a differential equation form. When we 
have disconnected vertices, we need write down all the possible configurations for the state of the particle at vertex 
[x, y) with probability of that possibility. For quantum percolation using D-DQW we find that the eight possible 
configuration can effectively describe the dynamics. We obtain the effective differential equation form by summing 
up the equations for each configuration with their respective probability of the possibility of occurrence. 



Stencil and the states for the vertex (x,y) 

In Fig. 7, we show the four configurations that result in transport of state from (x ± 1, y — 1) and (x, y — 1) to (x, y). 
In Fig. 8, we show the other four configurations that results in transport of state from [x ± 1, y — 1) to (x, y — 1) and 
the configuration that gets trapped at position {x,y — 1). This is equivalent configurations that results in transport 
of state from (x ± l,y) to {x,y) and the configuration that gets trapped at position (x,y). 

Figure. 7(a): Completely connected state and the expression are same as Eqs. (6) and (7) in the paper: 

i=i.y = COs(ff)Vi + i, v _i " i M^l-l,y-l (15) 
1pt,y = cos ( d )^l-l,y-l ~ * sin Wx+l,l/-l ' ( 16 ) 
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o- 



G- 



(a) (b) (c) (d) 

FIG. 7: (Color online) Four possible configuration when the transport can happen to vertex (x, y). Green, red and blue arrows representing 
the direction of | 4-), It) an d both the states, respectively. Transport from (a) (x + 1, y — 1) and (x — 1, y — 1) to [x, y), (b) (x — 1,1/ — 1) 
and (x, y — 1) to (x, y), (c) (x + 1, y — 1) and (x, y — 1) to (x, y) and (d) (x, y — 1) to (x, j/). 



These equations can be easily decoupled and written as, 



lit) 



/,4-(t) 



/,4-(t) 



(17) 



Subtracting 2 [1 + cos(#)] V'i,y from both sides of the preceding expression, we obtain a difference form which can be 
written as a second order differential wave equation 



^-cos(<?)|^ + 2[l-cos(0)] 



W,l(t) = o 

TX,y yj ' 



(18) 



When p is the percentage of connection in the lattice the probability for the above possibility is p 3 



Figure. 7(b) : Because of the missing edge from vertex (x, y — 1) to (x + 1, y — 1), the state arriving from (x, y — 2) 
will remain in (x, j/ — 1) to which the state from [x — 1, y — 1) join to proceed further to (x, y). For this configuration 
the state at (x, y) will be : 



After decoupling we get 



i>iy = cos(fl)^ il ,_ 1 - isin(^)[^-i, y -i + <„_i 



(19) 
(20) 

(21) 



Subtracting both sides of the preceding expression by [2 + cos(#)]?/>2 y + ipl+i y we obtain a difference form which can 
be written as 



o2 Q 

_ +cos (0)_ + [2-3cos(0)] 

The right hand side can be further simplified to obtain 

d 2 



Yx,y 



d_ 

dy 



ay oyox ox ay 



ii, v = o. 



(22) 



(23) 



The probability for the above possibility is p 2 (l — p). 

Figure. 7(c): For this configuration the state at (x,y) will be : 

^l, y = ^m(6) Vi, y _i +*l>t+i, v -i + cos ( d )$l, y -i- 



(24) 
(25) 
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(a) (b) (c) (d) 

FIG. 8: (Color online) During the transport the state from (a) (x — 1, y — 1) gets trapped at (x, y — 1) and (b) (x + 1, y — 1) gets trapped 
at (x, y — 1). (c) Due to the absence of edge from (x, y — 1) to (x, y) the states from (x ± 1, y) move back to (x, ±1, y). (d) absence of any 
transport or trapping. 



After decoupling we get 



(26) 



Subtracting both sides of the preceding expression by [2 + cos(8)]ipl y + ipl-i v we obtain a difference form which can 



be written as 



^-cos(0)^ + [2-3c OS (0)] 



The right hand side can be further simplified to obtain 



lb 1 



dy 



A " 7HF t 1 + cos ( )]|- - I" + t 1 " 3cob(0)] 
ay z oyax ox oy 



The probability for the above possibility is p 2 (l —p). 
Figure. 7(d): 

V£ v = cos(0)V£ s _i -isin^)^.! 







* =0. 



After decoupling we get 



ipl, v = -isia.{9)ipi )y _ 1 + cos(6»)V'I iy _i. 



Subtracting both the sides by 2ipx)y and writing the difference form as a differential equation, 



9 2 
dy 



^3 + 2 [1 - cos(0)] 



^ = o. 



The probability for the above possibility is p(l — p) 2 . 
Figure. 8(a): 

^i,y-l = C0S (°H'i,y-l -ism(d) ^-1,2,-1 + i>l,y-l 



(27) 



(28) 



(29) 
(30) 

(31) 
(32) 



(33) 
(34) 
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which is equivalent to 



i'l.y = -isin(e)ipi >y + cos(0) ipl_ ltV + ip$ >y 



After decoupling we get 



2[1 - coa(fl)]^ = [cos(0) _ !]^T_ 



(35) 
(36) 

(37) 



Subtracting both sides of the expression by [cos(0) — lJi/4 y we get a difference form that can be written in the 
differential equation for as, 



[1 - cos(0)] 

The probability for the above possibility is p(l — p) 2 . 
Figure. 8(b): 



d_ 

dx 



i>l,v = o 



(38) 



which is equivalent to 



After decoupling we get 



i>i, y -i = cos(0) + ipi+i, y -i ~ is H°)i>l,y-i 

^l,y-l = ~i Sm (#) IPty-1 + ^t+l.y-1 + COS Wt,y-l 



i>l, y = -* sin(0) Ui y + Vi+i, J + cob(0)iI>I v 



2[l-cos(0)]^ = [cos(0)-l]v4 +lir 



(39) 

(40) 
(41) 

(42) 

(43) 
(44) 

(45) 



Subtracting both sides of the expression by [cos(0) — y we get a difference form that can be written in the 
differential equation for as, 



[1 - cos(0)] 

The probability for the above possibility is p(l — p) 2 ■ 
Figure. 8(c): 



d_ 

dx 



(46) 



which is equivalent to 



i>t,y-i = cos(0)^+i,y_i - isin(0)v4-i, y -i 
4>t, y -i = cos(0)Vl_i iI ,_i - ism(0)V4 + i,i,-i • 

i>i, y = coa(9)ipi +1>y -iam{6)ipl_ hy 
i>t, y = cos(0)^l_i, y - isin(0)^ +liS/ • 



(47) 
(48) 

(49) 
(50) 



After decoupling we get 



2^=cos(0) 



/,t(40 , /T(40 

™x-l,y "t" rx+l,y 



(51) 
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Subtracting both sides of the preceding expression by 2 cos(9)ipt,y and by writing the difference form to the differential 
equation form we get: 



cos(e)^-2[cosW-l]]^ ) =0. 
The probability for the above possibility is p 2 (l —p). 
Figure. 8(d): 

^t,y-l = C0S ( S )^i >V -l -isMOWty-l 

i>t, y -i = -* sin Wi,y-i +™W4-i- 

which is equivalent to 

V>*,„ = cos (°)i>i, v -iBm(e)ii>l tV 

i't, y = -ism(6)ip^ y + cos{e)ipl t y. 
These expressions can be decoupled and written as: 

2[l-cos(^W)=0. 



(52) 



(53) 
(54) 

(55) 
(56) 

(57) 



The probability for the above possibility is (1 — p) 3 . 

Adding up the differential equations for all eight configuration as a product of the probability of its possibility we 
get : 



p 3 + 2p 2 {l-p)+p{l-pf 



dy 2 



p 3 -p 2 (l-p)] cos(#)|^ - 2p 2 (l -p)[cos(0)^ + ^ 



+[l-cos(0)] 2p 3 + 5p 2 {l-p) + 8p(l-p) 2 + 2(l-p) 3 + P 2 {1 -p)[l -3cos(0)] 



1>lv = 0. 



(58) 



Similar expression can be obtained for ijj^ . 



